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ABSTRACT 

We present an empirical model of Comptonization for fitting the spectra of 
X-ray binaries. This model, simpl, has been developed as a package imple- 
mented in XSPEC. With only two free parameters, simpl is competitive as the 
simplest empirical model of Compton scattering. Unlike other empirical mod- 
els, such as the standard power-law model, simpl incorporates the basic physics 
of Compton scattering of soft photons by energetic coronal electrons. Using a 
simulated spectrum, we demonstrate that simpl closely matches the behavior 
of physical Comptonization models which consider the effects of optical depth, 
coronal electron temperature, and geometry. We present fits to RXTE spectra 
of the black-hole transient HI 743-322 and a BeppoSAX spectrum of LMC X-3 
using both simpl and the standard power-law model. A comparison of the re- 
sults shows that simpl gives equally good fits and a comparable spectral index, 
while eliminating the troublesome divergence of the standard power-law model 
at low energies. Importantly, simpl is completely flexible and can be used self- 
consistently with any seed spectrum of photons. We show that simpl - unlike 
the standard power law - teamed up with diskbb (the standard model of disk 
accretion) gives results for the inner-disk radius that are unaffected by strong 
Comptonization, a result of great importance for the determination of black hole 
spin via the continuum-fitting method. 

Subject headings: Astrophysical Data: Data Analysis and Techniques 
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Introduction 



Spectra of X-ray binaries typically consist of a soft (often blackbody or bremsstrahlung) 
component and a higher-energy tail component of emission, which we refer to generically 
as a "power law" throughout this work. The origin of the power-law component in both 
neutron-star and black-hole sy stems is widely at t ributed to Compton up-scatter ing of soft 
photons by coronal electrons (IWhite et al.l Il995l ; iRemillard fe McClintockl 120061 . hereafter 
RM06). This component is present in the spectra of essentially all X-ray binaries, and it 
occurs for a wide range of physical conditions. 

The tail emission is generally modeled by adding a simple power-law component to 
the spectrum, e.g., via the model powerlaw in the widely used fitting package XSPEC 
(jArnaudlll996l ). A few of the many app lications where po wer-law models are employed in- 



clude: mo deling the thermal continuum (IShafee et al.ll2006l ) or the relativistically-broadened 



Fe K line ( [Miller et al.l 120081 ) in order to obtain estimates of black- hole spin; modeling the 



(White & Holt : 


L982) 


Gierlihski et al. 


1999) 



9821) or a substantial corona that scatters photons up to MeV energies 



19991 ); and classifying patterns of distinct X-ray states, e.g., in black- hole 



binaries (RM06). 

Because of the importance of the power- law component, several physical models have 
been developed to infer the conditions of the hot plasma that cause s the Comptoniz ation. 

EQ- 



Mode ls of this va riety that ar e available in XSPE C are compTT (ITitarc huk 



PAIR (CoppI 199S), cq mpTB ( Farinelli et al 



( INishimura et al 



198 6J), thcomp (IZycki et al 



2008), bmc (ITitarchuk et al.lll997l). co mpbb 



19991 ). compls (lLamb fe Sanfordlll979l ). and 



1994ft 



COMPPS (IPoutanen fe Svenssonlll996l ). It is essential to use such physical models when one 
is focused on understanding the physical conditions and structure of a scattering corona or 
other Comptonizing plasma. 

Often, however, the physical conditions of the Comptonizing medium are poorly un- 
derstood or are not of interest, and one is satisfied with an empirical model that seeks 
to match the data with no pretense that the model describes the physical system. The 
model powerlaw is one such empirical model which has been extraordi narily widely used 



in modeling black-hole and neutron-star binaries (see text & references in IWhite et al 



Tanaka fe Lewin 



2002 



1995 



Brenneman fe Revnoldsll2006l : RM06) and AGN (e.g 



1995; 



Zdziarski et al. 



Brenneman fe Reynolds! 2006 ). However, powerlaw introduces a serious flaw: at low 



energies it rises without limit. The divergence at low energies is unphysical, and it often 
significantly corrupts the parameters returned by the model component with which it is 
teamed (e.g., the widely used disk blackbody component diskbb; §3). 
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An excellent alternative to the standard power-law model for describing Compton scat- 
tering is provided by a convolution model that is based on a Green's function that was formu- 



lated decades ago (jShapiro. Lightman fc Eardleyll 19761 ; iRybicki fc Lightmanlll979l ; ISunyaev fc Titarchuk 
19801 ; iTitarchukl Il994l ). In this approach the power- law is generated self-consistently via 
Compton up-scattering of a seed photon distribution; consequently, the power-law naturally 
truncates itself as the seed distribution falls off at low energies. 

In this paper, we present our implementation of a flexible convolution model named 
SIMPL that can be used with any spectrum of seed photons. For a Planck distribution we show 
that simpl gives identical results to BMC, as expected since the two models are functionally 
equivalent (§2.3). Although simpl has only two free parameters, the same number as the 
standard powerlaw, this empirical model is nevertheless able to very successfully fit data 
simulated using compTT, a prevalent physical model of Comptonization (§2.2). 

We analyze data for two black hole binaries and illustrate the flexibility of simpl by 
convolvi ng simpl with piSKB B, the workhorse accretion disk model that has been used for 
decades ( iMitsuda et al.lll984l ). Our principal result is that simpl in tandem with diskbb 
enables one to obtain fitted values for the inner-disk radius R[ n for strongly- Comptonized 
data that are consistent with those obtained for weakly- Comptonized data (see §3.2). The 
standard power law, on the other hand, delivers very inconsistent values of _Ri n . As we 
show in §4.2, this result is consequential for the measurement of black hole spin via the 
continuum-fitting method: It implies that using simpl in place of the standard power law 
one can obtain reliable measurements of spin for a far wider body of data than previously 
thought possible, and for more sources (e.g., Cyg X-l). 

In §2 we outline the model and in §3 we present a case study with several examples. 
We discuss the implications and intended applications of the model in §4 and conclude with 
a summary in §5. 



2. The Model: simpl 



The model simpl (SIMple Power Law) functions as a convolution that converts a fraction 
of input seed photons into a power law (see eq. [I]). The model is currently available 
m XSPEC0. In addition to simpl- 2, whic h is ou r im plementation of the c l assica l model 
described by Ishapiro. Lightman fc Eardley ( 19761 ) and Sunyaev fc Titarchuk ( 1980 ). which 
corresponds to both up- and down-scattering of photons, we offer an alternative "bare-bones" 



http://heasarc.nasa.gov/xanadu/xspec/manual/XSmodelSimpl.html 
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implementation in which photons are only up-scattered in energy. The physical motivations 
behind the two versions of the model are described in §2.1, and the corresponding scattering 
kernels — the Green's functions — are given in equation (2) and equation (3), respectively. 

The parameters of simpl and the standard powerlaw model are similar. Their prin- 
cipal parameter, the photon index T, is identical. However, in the case of simpl the normal- 
ization factor is the scattered fraction /so, rather than the photon flux. The goal of simpl 
is to characterize the effects of Comptonization as simply and generally as possible. In this 
spirit, all details of the Comptonizing medium, such as its geometry (slab vs. sphere) or 
physical characteristics (optical depth, temperature ,thermal vs. non-thermal electrons ), 
which would require additional parameters for their description, are omitted. 

It is appropriate to employ simpl when the physical conditions of the Comptonizing 
medium are poorly understood or are not of interest. When the details of the Comptonizing 
medium are known, or are the main object of study, one should obviously use other models 
(e.g., compTT, compps, thcomp, etc.), which are designed specifically for such work. 
simpl, on the other hand, is meant for those situations in which a Compton power-law 
component is present in the spectral data and needs to be included in the model but is not 
the primary focus of interest, simpl should thus be viewed as a broad-brush model with the 
same utility as powerlaw but designed specifically for situations involving Comptonization. 

By virtue of being a convolution model, SIMPL mimics physical reprocessing by tying 
the power-law component directly to the energy distribution of the input photons. The 
most important feature of the model is that it produces a power-law tail at energies larger 
than the characteristic energy of the input photons, and that the power law does not extend 
to lower energies. This is precisely what one expects any Compton-scattering model to do 
and is a general feature of all the physical Comptonization models mentioned above. In 
contrast, the model powerlaw simply adds to the spectrum a pure power-law component 
that reaches all the way downward to arbitrarily low energies. The difference between simpl 
and powerlaw is thus most obvious at soft X-ray bands where simpl cuts off in a physical ly 



natural way whereas powerlaw continues to rise without limit (e.g., see lYao et al.ll2005l ). 



Two assumptions underlie simpl. The first is that all soft photons have the same 
probability of being scattered (e.g., the Comptonizing electrons are distributed spatially 
uniformly). This is a reasonable assumption when one considers that, even in the best 
of circumstances, almost nothing is known about the basic geometry of the corona. For 
example, usually the corona is variously and crudely depicted as a sphere, a slab, or a lamp 
post. The second assumption is that the scattering itself is energy independent. This is 
again reasonable given the soft thermal spectra of the seed photons that are observed for 
black-hole and neutron-star accretion disks with typical temperatures of ~ 1 keV and a few 
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keV, respectively. For example, in the extreme case of a 180° back-scatter off a stationary 
electron, a 3 keV seed photon suffers only a 1% loss of energy, and even a 10 keV photon 
loses only 4% of its initial energy. 

Figure 1 shows sample outputs from simpl when th e input soft photons are modeled 



by the multi-temperature disk blackbody model diskbb (IMitsuda et al.lll984l ). Results are 
shown for both SIMPL-2 and SIMPL- 1, our alternative version of SIMPL that includes only 
up-scattering of photons; the spectra are shown for T = 2.5 and a range of values of /so 
Note the power-law tails in the model spectra at energies above the peak of the soft thermal 
input and the absence of an equivalent power-law component at lower energies. This is 
the primary distinction between simpl and powerlaw. SIMPL-2 and simpl- 1 give similar 
spectra, but the spectrum from simpl- 1 has a somewhat stronger power-law tail for the 
same value of /so This is because simpl- 1 transfers all the scattered photons to the high 
energy tail, whereas simpl-2 has double-sided scattering. Therefore, for the same value of 
/sc, fewer photons are scattered into the high-energy tail with simpl-2. Correspondingly, 
when fitting the same data, simpl-2 returns a larger value of fsc compared to simpl- 1 (for 
examples, see §3 and Table 2). 



2.1. Green's Functions 



Given an input distribution of photons n- m (Eo)dEo as a function of photon energy E , 
simpl computes the output distribution n out {E)dE via the integral transform: 



He, 



:(E)dE 



hc)n m {E)dE + / 



sc 



n in (E )G(E;E )dE 



dE. 



X 



A fraction (1 — /sc) of the input photons remains unscattered (the first term on the right), 
and a fraction / sc is scattered (the second term). Here, E min and E max are the minimum 
and maximum photon energies present in the input distribution, and G(E; Eq) is the energy 
distribution of scattered photons for a 5-function input at energy E , i.e., G(E;E ) is the 
Green's function describing the scattering. 

We now describe the specific prescriptions we use for simpl-2 and simpl- 1. We also 
discuss the physical motivations b ehind these prescripti o ns, dr awing heavily on the theory 
of Comptonization as described by lRybicki &: Lightmanl (ll979l . hereafter RL79). 



- 6- 



2.1.1. SIMPL-2 



In sec. 7.7, RL79 discuss the case of unsaturated repeated sca ttering by nonrelativistic 
thermal electrons. Following IShapiro. Lightman fc Eardleyl (119761 ). they solve the Kompa- 
neets equation and show that Comptonization produces a power-law distribution of photon 
energies (eq. 7.76d in RL79). There are two solutions for the photon index V: 



r 2 



9 4 

4 + ? 

9 4 

4 + ? 



where the Compton y parameter is given by y = (4fcT e /m e c 2 )Max(r es , r c 2 s ). Up-scattered 
photons have a power-law energy distribution with photon index I\ and down-scattered 
photons have a different power-law distribution with photon index IV 

We model this case of nonrelativistic electrons with the following Green's function 
(jSunyaev fc Titarchuldll980l ; lTitarchuklll994! ; lEbisawalll999l ). which corresponds to the model 
SIMPL-2: 



G(E; E )dE 



(r-l)(r + 2) ((E/E )- r dE/E , 
(1 + 2T) \(E/E ) r+1 dE/E , 



E > Eq 
E < Eq. 



(2) 



The function is continuous at E = Eq, is normalized such that it conserves photons, and 
holds for all T > 1. Substituting <^fy in §B) we see that simpl-2 has two parameters: fsc and 
r. Although the model makes use of two power laws, their slopes are not independent. 

As in the case of the standard power law, simpl includes no high energy cutoff. Tech- 
nically, for any complete model of Comptonization, the up-scattered power-law distribution 
is cut off for photon energies larger than kT e . To avoid increasing the complexity of our 
model, we have ignored this detail; extra parameters could easily be added to account for 
high energy attenuation if desired. By keeping the model very basic, simpl is a direct two- 
parameter replacement for the standard power law while bridging the divide between the 
latter model and physical Comptonization models. 



2.1.2. SIMPL- 1 

The Green's function (2) is obtained by solving the Kompaneets equation, which assumes 
that the change in energy of a photon in a single scattering is small. This assumption is not 
valid when the Comptonizing electrons are relativistic. 
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In sec. 7.3 of their text, RL79 discuss Compton scattering by relativistic electrons with 
a power-law distribution of energy: n e (E e )dE e oc E~ p dE e . In the limit when the optical 
depth is low enough that we only need to consider single scattering, they show that the 
Comptonized spectral energy distribution (SED) is a power law of the form P(E)dE oc 
i^-O-i)/ 2 . Equivalently, the photon energy distribution takes the form n(E)dE oc E~ r , with 
a photon index T = (p + l)/2. Hardly any photons are down-scattered in energy. 

In sec. 7.5, RL79 show that repeated scatterings produce a po wer-law SED even when 



the r elativistic electrons have a non-power- law distribution (see also iTitarchuk fc Lyubarskij 



19951 ). In terms of the mean amplification of photon energy per scattering A and the optical 
depth to electron scattering r es , the Comptonized photon energy distribution takes the form 
n(E)dE oc E~ r with a photon index T — 1 — lnr es / In A. For the specific case of a thermal 
distribution of electrons with a relativistic temperature kT e m e c 2 , the amplification factor 
is given by A = \Q(kT e / m e c 2 ) 2 . Once again, hardly any photons are down-scattered. 

For both cases discussed above, Comptonization is dominated by up- scattering and 
produces a nearly one-sided power-law distribution of photon energies. This motivates the 
following Green's function, valid for Y > 1, which we refer to as the model SIMPL-1: 

f(r - l)(E/E )- r dE/E , E>E 
G(E;E )dE= < J\ / QJ / 0' _ o ^ 

[0, E < E . 

The normalization factor (r — 1) ensures that we conserve photons. 

Although simpl-1 is most relevant for relativistic Comptonization, it can also be used as 
a stripped-down version of SIMPL-2 for non-relativistic coronae. The reason is that the low- 
energy power-law (E/Eq) v+1 in equation (2) almost never has an important role. There is 
not much power in this component, and what little contribution it makes is indistinguishable 
from the input soft spectrum. Therefore, even for the case of nonrelativistic thermal Comp- 
tonization, for which the Green's function (2) is designed, there would be little difference if 
one were to use simpl-1 instead of simpl-2. 



2.2. Comparison to compTT 



To illustrate the performance of SI MPL relative to othe r Comptonization models, we have 
simulated a 2 x 10 6 -count BeppoSAX (IBoella et al.lll997l ) observation using the compTT 
model in XSPEC v!2.4.0x. 



For our source spectrum, we adopt disk geometry, a Wien distribution of seed photons 
at kT = 1 keV, and a hydrogen column density of iV H = 10 21 cm -2 . We set the optical 
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depth and temperature of the Comptonizing medium to r c = 2 and kT e = 40 keV. Our 
simulation uses the LECS, MECS, and PDS detectors on BeppoSAX, which span a wide 
energy range ~ 0.1 — 200 keV (for details on the instruments, see §3). The total number of 
counts in the simulated spectra (~2x 10 6 ) corresponds to a 3 ks observation of a 1 Crab 
source. 

We analyze the simulated data with a model consisting of a blackbody (bb) coupled with 
SIMPL. We refer to this model as SIMPL0BB (the <g> is to emphasize that simpl represents a 
convolution). The best fits achieved have reduced chi-squared values of xl — 1-00 (simpl- l) 
and xl = 1-06 (simpl-2). The fitted BB temperatures are respectively 1.14 ± 0.02 keV and 
1.29 ± 0.01 keV compared to 1 keV in the original compTT model. Figure 2 shows the fit 
using simpl- 1 and Table 1 lists the best-fit parameters for both models. 

In comparison, compbb, an alternative model of Compton scattering that assumes slab 
geometry, fits our simulated spectrum comparably well as simpl, with xl — 1-05 (Table 1). 
COMPBB returns the same temperature as simpl-2, kT hh = 1.29 ± 0.01 keV. Compared 
to the compTT progenitor, compbb gives similar estimates of the coronal temperature 
kT e and optical depth r c (Table 1). Even though compbb is a physically more realistic 
model of coronal scattering than SIMPL, it does not outperform SIMPL in terms of fitting the 
COMPTT-generated data. Meanwhile, the model bb+powerlaw performs quite poorly, 
yielding xl > 2. Parameters for this fit are given in Table 1. Note that the derived N-r using 
powerlaw is much higher than either the original value or those from fits with simpl. 

Though simpl is a purely empirical model, we see that it can deliver a remarkably 
successfully fit to data simulated using the physical model compTT. Even for a very cool 
corona with electron temperatures as low as kT e = 20 keV, which causes compTT to produce 
noticeable curvature in the high-energy spectrum, we find that simpl-2 and simpl- 1 achieve 
reasonable fits with xl < 1-2. 

A significant virtue of SIMPL relative to the physical Comptonization models in XSPEC 
is that simpl can be employed in conjunction with any source of seed photons. The physical 
models, on the other hand, are typically restricted to treating only one or two predefined 
photon distributions. One standard choice of continuum model that is widely used in fitting 
Comptonized accretion disks is diskbb+compTT. With simpl, one would instead employ 
the model simpl®diskbb. The latter not only generates the power law self-consistently via 
up-scattering of the seed photons, but it also has two fewer parameters. 
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2.3. Bulk Motion Comptonization 



The model BMC describes the Comptonization of blackbody seed photons by a converg- 
ing flow of isothermal gas t hat is freely falling toward a compact object, i.e., bulk motion 



Comptonization (see, e.g., IShrader fc Titarchukl Il998l ; iTitarchuk et al.l 119971 ). BMC is an 



alternative to coronal Comptonization models and is structured identically to simpl-2®bb; 
both models are specified with just four parameters. As a direct demonstration in XSPEC 
that siMPL-2(g)BB and BMC are identical, we analyzed our simulated BeppoSAX spectrum 
described above using both models. We found that the returned values of the column density 
Nh, the blackbody temperature kT, and the photon index T agreed in each case to four or 
more significant figures. 

BMC has been variously used to support claims that Compton scattering off in-falling gas 
within several gravitatio nal radii gives rise to the ob s erved high energy power la w in several 



black-hole binaries (e.g.. IShrader &: Titarchukl Il998l . Il999l ; iBorozdin et al.lll999l ). However 



this is only one interpretation of the model; simpl-2CS>bb, and therefore BMC, can equally 
be used to support a more standard model of coronal scattering (operating with uniform 
efficiency at all energies, see §2 and §2.1.2). Thus, although BMC is designed specifically to 
model relativistic accretion inflows, its function is actually quite general. 

One virtue of simpl relative to BMC is that simpl does depend upon discerning the 
nature of the Comptonizing region, be it corona, relativistic in-falling gas, or other. Another 
virtue of simpl is that it fully incorporates the utility of BMC while allowing complete flexi- 
bility in the choice of the spectrum of seed photons, e.g., simpl®diskbb is more appropriate 
for modeling Comptonization in accretion disks than BMC, which is hardwired to a Planck 
function. 

The theory of bu lk motion Comptonization is developed further and rigorously in 
Titarchuk et al.l (Il9971 ). This paper describes a Green's function that is more appropri- 
ate than the one used in BMC. A complete version of this Green's function is incorporated 
into the more sophisticated model compTB. However, this model is again limited to treating 
scattering from a predefined set of (blackbody-like) seed photon distributions and includes 
additional free parameters. We find that the fitting results obtained using this Green's func- 
tion are intermediate between those given by simpl- 1 and simpl- 2 so long as the temperature 
of the in-flowing electrons, T e , is above the observed energy range. 
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3. Data Analysis 



In this section, we apply simpl to a sample of observations to illustrate how simpl com- 
pares with powerlaw. To this end, we have selected two black-hole binaries, H1743-322 and 
LMC X-3. H1743-32 2 (hereafter H1743) is an especially pristine black-hole transient (see 
Remillard et al.ll2006l ) since, for much of its 2003 outburst, its spectrum can be satisfactorily 
modeled with just absorb ed (Ah ~ 2.2 x 10 22 cm -2 ) thermal-disk and power-law components 
(jMcClintock et al.ll2007bi hereafter M07). In particular, the 122 days of contiguous spec- 
tral data on which we focus do not require any additional components to accommodate the 
reflection or absorption features that are often present in the spectra of black hole binaries. 

The s pectra of H1 743 were acquired by the Rossi X-ray Timing Explorer (RXTE) PCU- 
2 module (ISwanklll999l ). RXTE's best-calibrated PCU detector, and were taken in "standard 
2" format. All spectra have been background subtracted and have typical exposure times 
~ 3000 s. The customary systematic error of 1% has been added to all energy channels. The 
resultant pulse-height spectra are analyzed from 2.8 — 25 keV using XSPEC vl2.4.0x (see 
M07 for further details). 

While RXTE provides good spectral coverage in hard X-rays (> 10 keV), which is 
most important for constraining the power-law component, it is not sensitive at low energies 
(< 2.5 keV). Therefore, RXTE data are generally insensitive to Ah- To complement the 
RXTE observations presented here, we have selected a BeppoSAX observation of LMC X-3, 
a persistent and pred ominantly thermal black-hole sou rce with a very low hydrogen column 
(AT H 



4 x 10 20 cm- 2 ; 



Page et al.ll2003l : lYao et ajj|2005j ) 



The BeppoSAX narrow-field instruments provide sensitive measurements spanning a 
wide range in energy, from tenths to hundreds of keV. The low-energy concentrator system 
(LECS) and the medium-energy concentrator system (MECS) probe soft fluxes, from ~ 
0.1—4 keV and ~ 1.5 — 10 keV, respectively. The phoswich detector system (PDS) is 
sensitive to hard X-rays from ~ 15 — 200 keV, and the high-pressure gas scintillation counter 
(HPGSPC) covers ~ 4 — 100 keV. In this analysis, we consider only the LECS, MECS, and 
PDS because the statistical quality of the HPGSPC data is relatively poor. 

In reducing BeppoSAX data, we have followed t he protocols given in the Cookbook for 
BeppoSAX NFI Spectral Analysis fjFiore et al.lll999l ). We use pipeline products and extract 
spectra from 8' apertures centered on LMC X-3 for both the LECS and (combined) MECS 
detectors. For the PDS, which is a simple collimated phoswich detector, we selected the 
fixed rise-time spectrum. In our analysis, we have used standard response matrices and 
included blank-field background spectra with the appropriate scalings. No pile-up correction 
is necessary. 
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3.1. Steep Power Law State 

About a third of the way through its nine- month outburst cycle, HI 743 repeatedly 
displayed spectra in the steep power-law (SPL) state that were devoid of absorption features. 
A salient feature of the SPL state is the presence of a strong power-law component of 
emission. (For a review of black-hole spectral states and a precise definition of the SPL state, 
see Table 2 and text in RM06.) Twenty-eight such featureless spectra were consecutively 
observed over a period of about three weeks (spectra #58-85; M07). We focus here on one 
representative spectrum, #77. In Figure 3 we show our fits and the associated unabsorbed 
models obtained using diskbb+powerlaw and siMPL(g>DlSKBB. Fitted spectral parameters 
are presented in Table 2. 

The quality of fit (as measured by xl) using either model is comparable. Nevertheless, 
there are distinct differences between the models. The fits with simpl have a ~ 50% larger 
disk normalization compared to powerlaw and a ~ 40% lower Ah (Table 2). The fit using 
powerlaw diverges at low energies, as revealed by removing photoabsorption from the 
fitted models (panels on the right in Fig. 3). The effect is quite severe and has no obvious 
physical explanation. In contrast, the fit using simpl is well behaved and the unabsorbed 
model is not divergent. 

3.2. Thermal Dominant State 

The key feature of the thermal dominant (TD) state is the presence of a totally dominant 
and soft (kT ~ 1 keV) blackbody-like component of emission that arises in the innermost 
region of the accretion disk. The TD state is defined by three criteria, the most relevant 
of which here is that the fraction of the total 2-20 keV unabsorbed flux in the thermal 
component is > 75%. For the full definition of this state, see Table 2 in RM06. 

Here we have chosen H1743 spectrum #91 which belongs to a sequence of ~50 featureless 
spectra (#86-136; M07) in the TD state. This spectrum has T ~ 2, which is somewhat 
harder than usual, but is otherwise typical of H1743's TD state. Spectral fit results are 
shown in Figure 4. In addition, in order to further illustrate for the TD state the differences 
between simpl and powerlaw at energies below the « 2.5 keV response cutoff of RXTE, 
we use a BeppoSAX observation of LMC X-3; our results are illustrated in Figure 5. This 
observation was carried out on 1996 November 28 with exposure times of 1.8, 4.5, and 2 ks 
respectively for the LECS, MECS and PDS. 

As in §3.1, we fit these data using diskbb+powerlaw and siMPL<g>DlSKBB. The 
best-fit spectral parameters are listed in Table 2. Due to a calibration offset between the 
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various BeppoSAX instruments, we follow standard procedure and fit for the normalization 
of the LECS and PDS relative to the MECS, the best-calibrated of the three. We adopt 
the canonical limits of 0.7 - 1 for LECS/MECS and 0.77 - 0.93 for PDS/MECS. These 
normalizations are included in the tabulated results. 

A comparison of the results obtained with powerlaw and simpl confirms the trends 
highlighted in §3.1, namely the differences in normalization and Ah- However, they are more 
modest here because the Compton component is weaker in the TD state. 

3.3. Comparison of simpl and powerlaw 

An examination of Table 2 reveals the following systematic differences in the derived 
spectral parameters returned when fitting with simpl vs. powerlaw: simpl yields (i) a 
stronger and softer thermal disk component, i.e., a larger normalization and lower kT*; (ii) 
a generally steeper power law component (larger T); and (iii) a systematically lower Ah- As 
we now show, all of these effects can be simply understood. 

Because powerlaw produces higher fluxes than simpl at low energies, it tends to 
suppress the flux available to the (soft) thermal component, namely diskbb in the examples 
given here. This explains why powerlaw tends to harden the diskbb component and 
to steal flux from it (i.e., reduce its normalization constant). Meanwhile, at low energies 
the powerlaw component predicts artificially high fluxes that, in order to conform to the 
observed spectrum, depress the value of T. These differences between simpl and powerlaw 
are most pronounced when the power law is relatively steep, i.e., typically when T > 3. 

Modest and reasonable values of Ah are returned in fits using simpl, as well as compTT 
and other Comptonization models, because the Compton tail is produced by the up-scattering 
of seed photons and there is no power- law component at low energies. In contrast, POW- 
ERLAW continues to rise at low energies, which forces A H to increase in order to allow the 
model to fit the observed spectrum. This systematic difference is apparent in our fit results 
for the HI 743 spectra and is especially prominent in the case of the LMC X-3 spectrum 
for which Ah differs by a factor of two. For HI 743, the discrepancy in Ah is much less for 
the TD spectrum than for the SPL spectrum because the SPL state has both a steeper and 
relatively stronger power-law component. 

We turn now to consider the DISKBB normalization constant, which is proportional to 
Rf n , the square of the inner disk radius (see footnotes to Table 2). For the pair of HI 743 
spectra, we see that the disk normalization obtained with powerlaw is ~35% smaller in 
the SPL state than in the TD state (Table 2), indicating that R in is smaller for the SPL 
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state. With simpl, on the other hand, there is no significant change in the normalization, 
and hence both the SPL and TD states can be modeled with a disk that has the same inner 
radius. The radius is constant because simpl recovers the original (unscattered) flux emitted 
by the disk, which POWERLAW cannot do. 

This crucial ability to unify the inner regions of the accretion disk in thermally active 
states exhibiting both high and low levels of Comptonization (i.e., TD and SPL states) paves 
the way for a full genera l relativistic analysis which c an formally link R in to black-hole spin 
(see discussion in §4.2p . iKubota fc Makishimal (120041 ) similarly identified a constant radius 
for the black hole binary XTE J1550-564 between the TD and SPL states in an analysis using 
the model diskbb + thcomp. Because thcomp is implemented as an additive (i.e., non- 
convolution) model, Kubota & Makishima had to employ an awkward and ad hoc procedure 
to obtain their result (see their Appen dix). Their work i mproved upon a similar result 
obtained for black-hole GRO J1655-40 (IKubota et al.ll200ll ). With simpl, the modeling is 
significantly easier. 



4. Discussion 

4.1. Black Hole X-ray States 

A standard method of classifying X-ray states in black hole binaries involves spectral 
decomposition into two primary components - a multi-temperature blackbody disk, diskbb, 
and a Compton power law, powerlaw (RM06). This method is compromised by the use of 
the standard power law when the photon index is large (r > 3). In this case, at low energies 
the flux from the power law can rival or exceed the thermal component and thereby pollute 
it. As discussed in §2, intrusion of the power-law component at low energies is fundamentally 
inconsistent with Compton scattering. 

This difficulty in classifying states, which is caused by the use of powerlaw, is remedied 
by the use of simpl because the latter model naturally truncates the power-law component at 
low energies. It is useful to consider the intrinsic differences between the two models and how 
they influence the classification of black-hole X-ray states. Using powerlaw, the thermal 
disk and tandem power-law emission are modeled independently. On the other hand, under 
simpl all photons originate in the accretion disk. Some of these disk photons scatter into a 
power law en route from the disk to the observer. As described in §3.3, fits employing simpl 
imply stronger disk emission and weaker Compton emission than those using powerlaw. 
As a result, state selection criteria would need to be modified for classification using simpl. 
This topic is beyond the scope of the present paper. 
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4.2. Application to the Measurement of Black Hole Spin 



During the past few decades, the masses of 22 stellar black holes have been measured, 17 
of which are found in transient black ho le binaries (RM06 ) . Recently, we have measured the 
spins of four of these stellar black holes ( Shafee et al. 2006[ McClintock et al. 2006 ; Liu et al. 



20081 ) by fi tting their continuum spectra to our fu lly relativistic model of an accretion disk 
KERRBB2 (ILi et al.ll2005t iMcClintock et aJj [2000 ) plus the standard power-law component 



POWERLAW. In the continuum- fitting method, sp in is measured by estimating the inner 
radius of the accretion disk R in ( jZhang et al.lll997l ). We identify R m with the radius of the 



innermost stable circular orbit -Risco, which is predicted by general relativity. 

To date, we have conservatively selected only TD data for analysis (§3.2). Meanwhile, 
the transient black hole binaries spend only a modest fraction of their outburst cycle in 
the TD state and are often found in the SPL or some intermediate state (see Figs. 4-9 in 
RM06). Since for each source we seek to obtain as many independent measurements of the 
spin parameter as possible, our sole reliance on TD data has been a significant limitation. 
Using simpl, the highly Comptonized SPL state is now able to provide estimates of spin, a 
matter to be discussed more fully in (J. Steiner et al. 2009, in preparation). Not only will 
this allow us to substantially increase the size of our data sample for many sources, it also 
will likely allow us t o obtain spin measurements fo r sources such as Cygnus X-l that never 
enter the TD state JMcClintock & Remillanjfaod ). 



Our reason for developing simpl was to improve our methods for analyzing TD-state 
data in order to determine more reliable values of black hole spin. That simpl now allows us 
to determine spins for SPL-state data was a serendipitous discovery and a major bonus. We 
were motivated to develop simpl because we have been hampered by the use of powerlaw 
in two ways. Fir st, in all of our work we have exclusively used TD spectral data (e.g., 
Shafee et al.ll2006l ). which is maximally free of the uncertain effects of Comptonization. The 
selection of these data is problematic because, as indicated in §4.1 above, it can be affected 
in unknown ways when the spectral index of the Compton component is large. 

Secondly, of greater concern is the potential adulteration in the TD state of the thermal 
component by the power-law component, which can have an uncertain and sizable effect 
on the fitted value of the black-hole spin parameter. In the case of a number of spectra 
with steep power-law components, we found that the fitted values of the spin parameter 
were affect ed by the contribution o f the power law flux at energies below ~ 5 keV; e.g., 
see § 4.2 in IMcClintock et al.l (120061 ). In order to mitigate this problem, we applied in turn 
two alternative models: compTT, and a standard power-law component curtailed by an 
exponential low-energy cutoff, expabsx powerlaw. The former model was unsatisfactory 
because we were unable to fit for reasonable values of both the coronal temperature kT e and 
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the optical depth r c . The latter model was likewise unsatisfactory because it requires the 
use of an arbitrary cutoff energy E c . 

We developed simpl in order to sidestep these difficulties and uncertainties. At low 
energies, the model truncates the power law in the same physical manner as compTT 
and other sophisticated Comptonization models, simpl self-consistently ties the emergent 
power-law flux to the seed photons in order to deliver the power-law component via coronal 
reprocessing. An application of simpl to the measurement of the spin of the black-hole 
primary in LMC X-l is described in L. Gou et al. (2009). 

5. Summary 

We present a new prescription for treating Comptonization in X-ray binaries. While no 
new physics has been introduced by this model, its virtues lie in its simplicity and natural 
application to a wide range of neutron-star and black-hole X-ray spectra, simpl offers a 
generic and empirical approach to fitting Comptonized spectra using the minimum number 
of parameters possible (a normalization and a slope), and it is valid for a broad range of 
geometric configurations (e.g., uniform slab and spherical geometries). The scattering of a 
seed spectrum occurs via convolution, which self-consistently mimics physical reprocessing of 
photons from, e.g., an accretion disk. In addition to this physically motivated underpinning, 
simpl remains as unassuming as powerlaw but without its troublesome divergence at low 
energies. 

Our model is valid for all T > 1. We have shown that simpl is able to provide a good fit 
to a demanding simulated data set, which was generated with the widely-used Comptoniza- 
tion model compTT. Furthermore, we have demonstrated that simpl and powerlaw give 
very comparable results when fitting spectral data in terms of quality-of-fit and spectral 
index (see Table [2]). This quality of performance holds true not only for spectra with weak 
Compton tails (TD state) but also for spectra requiring a large Compton component (SPL 
state). In the latter case, the model based on simpl gives physically more reasonable results 
for the soft end of the spectrum (e.g., see §3.3). 

Using simpl®diskbb it will be important to revisit the classification of black hole states 
(RM06) for two reasons. First, the selection of TD data will no longer be adversely affected 
by the presence of a steep power-law component. Secondly, this model will allow some degree 
of unification of the TD state and SPL state, the latter being a more strongly Comptonized 
version of the former. In determining black hole spin via the continuum-fitting method using 
KERRBB2, simpl is a significant advance on three fronts: It will (1) enable the selection of 
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data with a dominant thermal component that is not mucked up by the effects of a divergent 
power-law component; (2) allow reliable spin measurements to be obtained using strongly- 
Comptonized SPL data, thereby substantially increasing the data sample for a given source; 
and (3) likely make possible the determination of the spins of some black holes that do not 
enter the TD state (e.g., Cygnus X-l). 

The authors would like to thank George Rybicki for discussions on the physics of Comp- 
tonization as well as Jifeng Liu, Lijun Gou, Rebecca Shafee, and Ron Remillard for their 
input on SIMPL. JFS thanks Joey Neilsen for enthusiastic discussions as well as comments 
on the manuscript, Ryan Hickox for suggestions which improved this paper, and Keith Ar- 
naud for helping implement simpl in XSPEC. The authors thank Tim Oosterbroek for his 
indefatigable assistance with the BeppoSAX reduction software. JFS was supported by the 
Smithsonian Institution Endowment Funds and RN acknowledges support from NASA grant 
NNX08AH32G and NSF grant AST-0805832. JEM acknowledges support from NASA grant 
NNX08AJ55G. 



A. XSPEC Implementation 

SIMPL is presently implemented in XSPEC. This version includes three parameters (two 
that can be fitted), the power-law photon index (r), the scattered fraction (/sc), and a 
switch to set up-scattering only (simpl- 1: switch > 0) and double-sided scattering (simpl- 
2: switch < 0). Since SIMPL redistributes input photons to higher (and lower) energies, for 
detectors with limited response matrices (at high or low energies), or poor resolution, the 
sampled energies should be extended or resampled within XSPEC to adequately cover the 
relevant range. For example, when treating the RXTE data in §3, which has no response 
defined below 1.5 keV, the command "energies 0.05 50 1000 log" was used to explicitly extend 
and compute the model over 1000 logarithmically spaced energy bins from 0.05 — 50 keV. 

Using simpl can be problematic when V is large, especially if the power-law component 
is faint or the detector response extends only to ~ 10 keV (e.g., Chandra, XMM or ASCA). 
When the photon index becomes sufficiently large, a runaway process can occur in which Y 
steepens and the scattered fraction becomes abnormally high (typically > 50%, inconsistent 
with a weak power law). This occurs because scattering redirects photons from essentially a 
5-function into a new function with characteristic width set by V. If V reaches large values 
(> 5), the scattering kernel will also act like a 5-function, and the convolved spectrum will 
be nearly identical to the seed spectrum. 
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In such circumstances, we recommend bracketing T. In practice, the power-law spectral 



indic es of black- hole systems are found to lie in the range 1.4 < T < 4 (IRemillard fc McClintock 



20061 ). An upper limit of T ~ 4 — 4.5 is typically sufficient to prevent this runaway effect, 



and this constraint should be applied if it is deemed appropriate for the source in question. 

We advise against applying simpl to sharp components such as spectral lines or reflec- 
tion components. SIMPL broadens spectral lines, and for prominent higher energy emission 
features (such as Fe Ka), simpl can saturate the power- law flux with scattering from the line 
itself. To prevent this from occurring, such components should be invoked outside the scope 
of simpl. For example, the XSPEC model declaration "model phabsx(smedgexsimpl(kerrbb)+laor)' 
would satisfy this recommendation. 



REFERENCES 

Arnaud, K. A. 1996, in ASP Conf. Series, Vol. 101, Astronomical Data Analysis Software 
and Systems V, ed. G. H. Jacoby & J. Barnes (San Francisco: ASP), 17 

Boella, G., Butler, R. C, Perola, G. C, Piro, L., Scarsi, L., & Bleeker, J. A. M. 1997, A&AS, 
122, 299 

Borozdin, K., Revnivtsev, M., Trudolyubov, S., Shrader, C, & Titarchuk, L. 1999, ApJ, 
517, 367 

Brenneman, L. W., & Reynolds, C. S. 2006, ApJ, 652, 1028 

Cowley, A. P., Crampton, D., Hutchings, J. B., Remillard, R., & Penfold, J. E. 1983, ApJ, 
272, 118 

Coppi, P. S. 1999, High Energy Processes in Accreting Black Holes, ed. J. Poutanen & 
R. Svensson (San Francisco: ASP), 161, 375 

Ebisawa, K. 1999, in ASP Conf. Series, Vol. 161, High Energy Processes in Accreting Black 
Holes, ed. J. Poutanen & R. Svensson (San Francisco: ASP), 39 

Farinelli, R., Titarchuk, L., Paizis, A., & Frontera, F. 2008, ApJ, 680, 602 

Fiore, F., Guainazzi, M., & Grandi, P. 1999, Cookbook for BeppoSAX NFI Spectral Analysis, 
Bepposax SDC 

Gierlihski, M., Zdziarski, A. A., Poutanen, J., Coppi, P. S., Ebisawa, K., & Johnson, W. N. 
1999, MNRAS, 309, 496 



-18- 



Gou, L., McClintock, J. E., Liu, J., Narayan, R., Steiner, J. F., Remillard, R. A., Orosz, 
J. A., k Davis, S. W. 2009, arXiv:0901.0920 [astro-ph] 

Kubota, A., Makishima, K., & Ebisawa, K. 2001, ApJ, 560, L147 

Kubota, A., & Makishima, K. 2004, ApJ, 601, 428 

Lamb, P. & Sanford, P. W. 1979, MNRAS, 188, 555 

Li, L.-X., Zimmerman, E. R., Narayan, R., & McClintock, J. E. 2005, ApJS, 157, 335 

Liu, J., McClintock, J. E., Narayan, R., Davis, S. W., & Orosz, J. A. 2008, ApJ, 679, L37 

McClintock, J. E., & Remillard, R. A. 2006, in Compact Stellar X-ray Sources, ed. W. Lewin 
& M. van der Klis (Cambridge: Cambridge Univ. Press), 157 

McClintock, J. E., Remillard, R. A., Rupen, M. P., Torres, M. A. P., Steeghs, D., Levine, 
A. M., & Orosz, J. A. 2007b, arXiv:0705.1034vl [astro-ph] 

McClintock, J. E., Shafee, R., Narayan, R., Remillard, R. A., Davis, S. W., & Li, L.-X. 2006, 
ApJ, 652, 518 

Miller, J. M., Reynolds, C. S., Fabian, A. C, Cackett, E. M., Miniutti, C, Raymond, J., 
Steeghs, D., Reis, R., & Homan, J. 2008, ApJ, 679, LI 13 

Mitsuda, K., Inoue, H., Koyama, K., Makishima, K., Matsuoka, M., Ogawara, Y., Suzuki, 
K., Tanaka, Y., Shibazaki, N., & Hirano, T. 1984, PASJ, 36, 741 

Nishimura, J., Mitsuda, K., & Itoh, M. 1986, PASJ, 38, 819 

Orosz, J. A. 2003, in IAU Symposium, Vol. 212, A Massive Star Odyssey: From Main Se- 
quence to Supernova, ed. K. van der Hucht, A. Herrero, & C. Esteban (San Francisco: 
ASP), 365 

Page, M. J., Soria, R., Wu, K., Mason, K. O., Cordova, F. A., & Priedhorsky, W. C. 2003, 
MNRAS, 345, 639 

Poutanen, J. & Svensson, R. 1996, ApJ, 470, 249 

Remillard, R. A. & McClintock, J. E. 2006, ARA&A, 44, 49 

Remillard, R. A., McClintock, J. E., Orosz, J. A., & Levine, A. M. 2006, ApJ, 637, 1002 

Rybicki, G. B. & Lightman, A. P. 1979, Radiative processes in astrophysics (New York, 
Wiley-Interscience) 



-19- 



Shafee, R., McClintock, J. E., Narayan, R., Davis, S. W., Li, L.-X., & Remillard, R. A. 2006, 
ApJ, 636, L113 

Shapiro, S. L., Lightman, A. P., & Eardley, D. M. 1976, ApJ, 204, 187 

Shrader, C. & Titarchuk, L. 1998, ApJ, 499, L31 

Shrader, C. R. & Titarchuk, L. 1999, ApJ, 521, L121 

Swank, J. H. 1999, Nuclear Physics B Proceedings Supplements, 69, 12 

Sunyaev, R. A. & Titarchuk, L. G. 1980, A&A, 86, 121 

Tanaka, Y. & Lewin, W. H. G. 1995, in X-ray Binaries, ed. W. Lewin, J. van Pardijs, & E. 
van den Heuvel (Cambridge: Cambridge Univ. Press), 126 

Titarchuk, L. 1994, ApJ, 434, 570 

Titarchuk, L. & Lyubarskij, Y. 1995, ApJ, 450, 876 

Titarchuk, L., Mastichiadis, A., & Kylafis, N. D. 1997, ApJ, 487, 834 

White, N. E. & Holt, S. S. 1982, ApJ, 257, 318 

White, N. E., Nagase, F., & Parmar, A. N. 1995, in X-ray Binaries, ed. W. Lewin, J. van 
Pardijs, & E. van den Heuvel (Cambridge: Cambridge Univ. Press), 1 

Yao, Y., Wang, Q. D., & Nan Zhang, S. 2005, MNRAS, 362, 229 

Zdziarski, A. A., Leighly, K. M., Matsuoka, M., Cappi, M., & Mihara, T. 2002, ApJ, 573, 
505 

Zhang, S. N., Cui, W., & Chen, W. 1997, ApJ, 482, L155 
Zycki, P. T., Done, C, & Smith, D. A. 1999, MNRAS, 309, 561 



This preprint was prepared with the AAS IATgX macros v5.2. 



Table 1. Results of Fitting a Simulated compTT Spectrum 



MODEL 




(10 22 cm- 2 ) 


r 


he 


Norm(PL) a 


kT 
(kcV) 


Norm b 


kT e 
(keV) 


Tc 


compTT c 




0.1 








1. 


0.001 


40. 


2. 


SIMPL-1 <g> BB 
SIMPL-2 ® BB 
COMPBB 
BB+POWERLAW 


1.00/731 
1.06/731 
1.05/731 
2.02/731 


0.28 ±0.01 
0.31 ±0.01 
0.31 ±0.01 
0.68 ±0.01 


1.41 ±0.02 
1.37 ±0.02 

1.00 ±0.01 


0.84 ±0.01 
0.87 ±0.01 


(5.0 ±0.2) x 10~ 3 


1.142 ± 0.015 
1.292 ± 0.010 
1.292 ± 0.010 
1.700 ± 0.008 


10.9 ±0.4 
7.8 ± 0.3 
19.7 ±0.7 
0.89 ± 0.02 


43.6 ±2.2 


2.21 ± 0.03 



a POWERLAW normalization given at 1 keV in photons s 'cm 2 keV 1 . 

b BB and COMPBB normalization = ( rj^it^kpc ) ^ or a blackbody of radius R at a distance D; COMPTT normalization is undefined. 
c COMPTT model set to disk geometry (geometry switch = 1). 



Table 2. Spectral Fit Results 



POWER LAW 



Source 
State 



Mission 
Detector 



MJD 



^disk a 
^Edd 



Ver.' 



he 



(10 2 



cm 



2 ) 



fcT* 
(keV) 



Norm c 



Norm(PL) d 



H1743 RXTE 52797.6 0.64/44 0.19 1.94 ±0.17 ■•■ ■■■ ■■• 1.189 ±0.011 568 ± 31 
SPL PCA 0.64/44 0.26 1.16 ±0.17 SI 2.65 ± 0.02 0.170 ±0.004 1.157 ±0.011 875 ± 49 
0.64/44 0.27 1.19 ±0.17 S2 2.65 ± 0.02 0.224 ± 0.005 1.162 ±0.011 878 ± 50 

H1743 RXTE 52811.5 0.67/44 0.22 1.53 ±0.15 ■•■ ■■• ■■• 1.106 ±0.006 869 ± 33 

TD PCA 0.67/44 0.23 1.48 ±0.15 SI 1.98 ± 0.03 0.030 ± 0.001 1.104 ±0.006 909 ± 35 

0.67/44 0.23 1.48 ±0.15 S2 1.98 ± 0.03 0.037 ± 0.001 1.105 ±0.006 910 ± 35 



2.64 ±0.02 9.92 ±0.54 



1.98 ±0.03 0.52 ±0.05 



LMCX-3 BeppoSAX* 50415.5 1.05/729 0.58 0.073 ± 0.008 ■•■ ■■• ■■• 1.279 ±0.011 24.5 ± 0.8 2.19 ±0.11 0.055 ± 0.010 

TD LECS.MECS, 1.08/729 0.60 0.044 ± 0.003 SI 2.41 ± 0.45 0.062 ± 0.021 1.238 ±0.013 30.4 ±1.2 

PDS 1.08/729 0.59 0.044 ± 0.003 S2 2.46 ± 0.48 0.085 ± 0.033 1.239 ±0.012 30.3 ±1.1 ■■• 



a Bolometric (0.1 — 20 keV) luminosity of t he disk component i n Eddington units. For H1743, we adopt nominal values: M = 10 Mq, D = 7.5 kpc, and i = 70°. The fiducial values used 
for LMC X-3 are M = 7.5 M and % = 67° l lCowlev et al.lll983tlOrosdl2003t) . For fits using SIMPL, this quantity describes the seed spectral luminosity. 

b Version of SIMPL being used, i.e., SI for SIMPL-1 and S2 for SIMPL-2. 

c For an accretion disk inclined by i to the line of sight, with inner radius ilj n at distance D, Norm = ( d/io k*pc ) cos *- 
d POWERLAW normalization given at 1 keV in photons s _1 cm~ 2 keV - 1 . 

G The cross-normalizations for C LM = LECS/MECS and C PM = PDS/MECS are fitted from 0.7 - 1 and 0.77 - 0.93 respectively. C LM = 0.802 ± 0.283, 0.814 ± 0.008, 0.813 ± 0.008 for 
the fits with POWERLAW, SIMPL-1, and SIMPL-2. Cpm is pegged at 0.93 for the same fits. 



Note. - 



All errors are presumed Gaussian and quoted at la. 
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Fig. 1. — Spectral energy density vs. photon energy for a sample spectrum calculated with 
SlMPL-1 (solid lines) and simpl-2 (dashed lines). The models conserve photons and Comp- 
tonize a seed spectrum, which in the case shown is diskbb with kT* = 1 keV (black line). 
Ascending colored lines show increasing levels of scattering, from /sc = 1 — 100%. 
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Fig. 2. — The data correspond to a simulated BeppoSAX observation with a total of 2.1 xlO 6 
counts; the spectrum was generated using compTT. The histogram shows the fit achieved 
using simpl-1. This fit is performed over the recommended energy ranges of the narrow-field 
instruments (NFI), as given by the Cookbook for BeppoSAX NFI Spectral Analysis, yielding 
xt = 1.00. For details, see Table 1. This example demonstrates the ability of simpl to 
match a representative spectrum generated by a physical model of Comptonization. 
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Fig. 3. — left: Unfolded spectral fits to an RXTE observation of H1743 in the SPL 
state and right: the corresponding unabsorbed models. Data are fitted using (a,b): 
PHABSx(DISKBB+POWERLAW), (c,d): PHABSx(siMPL-l(g)DISKBB), (e,f): PHABS X (SIMPL- 
2®diskbb). The composite model is represented by a solid black line and the emergent disk 
and Compton components are shown as red and blue dashed lines respectively. The seed 
spectrum for simpl is shown (dashed) in green. Contrasting behaviors between simpl and 
powerlaw are most clearly revealed in the unabsorbed models at low energies. Spectral 
parameters are given in Table 2. 



-25 - 




1 2 3 4 5 10 20 3040 1 2 3 4 5 10 20 30 40 

Energy (keV) 



Fig. 4. — Same as Figure 3 except that the results shown here are for an RXTE observation 
of HI 743 in the TD state. The systematic differences between the simpl and powerlaw 
fits are greatly reduced compared to the differences shown for the SPL example in Figure 3. 
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Fig. 5. — Same as Figures 3 and 4 for a TD BeppoSAX spectrum of LMC X-3. The data 
have been rebinned for plotting purposes only and both LECS and PDS counts have been 
rescaled by the fitted normalizations given in Table 2. At low energies (below ~0.5 keV), 
the unabsorbed model is strongly compromised for fits with powerlaw 



